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ABSTRACT 


The thermal electron heating rate, Q e , is an important heat source term 

in the ionospheric electron energy balance equation, representing heating by 

photoelectrons or by precipitating higher energy electrons. A formula for the 

thermal electron heating rate is derived from the kinetic equation using the 

electron-electron collision operator as given by the unified theory of Kihara 

and Aono. This collision operator includes collective interactions to produce 

a finite collision operator with an exact Coulomb logarithm term. The derived 

heating rate Q e is the sum of three terms, Q e = Q p + S -f Q int , which are 

respectively: 1) primary electron production term giving the heating from 

newly created electrons that have not yet suffered collisions with the ambient 

electrons, 2) a heating term evaluated on the energy surface m v^2 * E_ at the 

6 1 

transition between Maxwellian and tail electrons at E^, and 3) the integral 
term representing heating of Maxwellian electrons by energetic tail electrons 
at all energies > E_ Published ionospheric electron temperature studies have 

i • 

used only the integral term Q int with differing lower Integration limits. 

There can be a significant numerical difference between Q e and Q^ nt . Use of 
the imcomplete heating rate could lead to erroneous conclusions regard' 
electron heat balance, since Q e is greater than t by as much as a factor of 
two. The sensitivity of the heating rate to the method of calculating the 
energetic (tail) electron distribution function, using either a linear or a 
quadratic collision operator is demonstrated. Choice of the transition 
energy, E-j., between the thermal (Maxwellian) population and the energetic tail 
electrons significantly affects the magnitude of the individual heating rate 
terms. The net heating rate Q e is less sensitive to the value of E T then are 


the individual terms 


1 . INTRODUCTION 


M 


'The Earth's ionospheric layer has been studied ever since its importance 
for long distance radio communication was discovered. The earliest iono- 
spheric studies dealt mainly with its density variations with altitude, time 
and location. The thermal properties of the Earth's ionosphere first came 
under Intensive study in the late 1950' s and early 1960's when rockets became 
available for in situ measurements. The first in situ electron temperature 
measurements were made by Boggess et al. (1959), Spencer et al. (1962), and 
Brace et al. (1963); these measurements clearly demonstrated that the electron 
temperatures exceeded the neutral temperature in the dayside ionosphere. 
Predating these measurements by about 15 years, Drukarev (1946) predicted the 
existence of hot ambient electrons in the ionosphere due to heating from 
energetic photoelectrons. Early mod°l calculations of ionospheric electron 
temperatures were made by Hanson and Johnson (1961), Dalgarno et al. (1963), 
and Hanson (1963). The$?e calculations were made by balancing the ambient 
electron heating rate with the rate of electron cooling from collisions with 
ions and neutral particles. In their early work, Hanson and Johnson (1961) 
evaluated the ambient electron heating rate as the product of a constant 
heating efficiency multiplied by the total photoelectron production rate. 

Electron heat balance studies after this pioneering work have calculated 
the heating rate as the product of the energetic electron (tail) distribution 
function and the energy loss rate between tail electrons and Maxwellian 
electrons, integrated over the energy range of the tail electrons. This 
formula for calculating the heating rate is an assumed formula. The deriva- 
tion of the heating rate presented in this paper will show that it is one of 
the three terms contributing to the heating rate. The other two terms are 
significant and therefore can not be neglected. 
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The reader is referred ';o the review paper by Schunk and Nagy (1978) for 
a general discussion of ionospheric electron temperature theory and observa- 
tions, and for further references. The remainder of this paper is devoted to 
the thermal electron heating rate, how it is derived, a formula for it, and 
some numerical examples. In section 2 we give the basic equations (the 
electron kinetic equation and the electron fluid equation) from which the 
ambient heating rate is derived. An alternative electron heat balance 
equation is discussed and the need for partitioning the electrons into two 
populations is given. The formula for the ambient electron heating rate is 
derived in section 3 and some numerical examples and comparisons with earlier 
heating rates are given in section A. Concluding remarks are given in section 
5. Appendix A sketches a derivation of the electron-electron collision 
operator which is required in the heating rate derivation. 


2. BASIC EQUATIONS 

The starting point in the derivation of the thermal electron heating rate 
is the kinetic equation for the electron velocity distribution function, f(v). 
We write the kinetic equation in the form (Burgers, 1969): 

_ +v . + .? mf = p + _ (i) 

3x 3 v 

where $ is the force per mass, p is the electron production rate from photo- 
ionization of neutrals, and 6f/6t is the collisional rate of change of f. 

The electron heating rate is a term in the electron energy balance 
equation, which is obtained by taking the energy moment of the electron 
kinetic equation. If we follow the usual procedure of integrating over all 
velocities (Including both populations of thermal and tail electrons) we 
obtain the energy equation: 
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( 2 ) 




»eW 
+ div q„ 


3/2 kT e + N e kT e dtv E e + E e : — £ 
- <(E - | kT e ) p > + <(E - | kT e ) 



2 

mu 

where < > denotes the Integral over all velocities v, E ■ — is the thermal 

energy, u « v - i$ , N e is the electron density, T e is the electron tempera- 

ture, t is the average electron velocity, ? is the traceless momentum 

tensor, and q is the electron heat flux, 
e 

This equation is consistent with the energy equations of Shkarofsky et 
al. (1966) and Burgers (1969). The form appears different from the latter two 
references because the density, N e , has been removed from the substantial 
derivative and we include the electron production term p. The density, 
temperature, average velocity, momentum tensor, and heat flux are defined (in 
Eq. 2) for the total electron population, since the moments of £ were taken 
over the entire velocity space (see Chapman and Cowling, 1960 for the 
definitions) . 

The first term on the right side of the energy balance equation repre- 
sents direct heating from electrons which are created by ionization of neutral 
particles. The second term is a net cooling rate from the totality of 
electron-ion and electron-neutral collisions. There is no thermal electron 
heating rate term in this equation because electron-electron collisions 
conserve energy and particles, and the integrals were carried out over all 
velocities. The heating of the thermal electron population (the Maxwellian 
electrons) balances the cooling of the energetic electron population (the 
enhanced tail electrons) so that the sum is zero. Thus the usual hydrodynamic 
energy equation, Eq. 2, does not contain the heating rate term that we wish to 
derive. Also, the transport terms in Eq. 2 such as div q g contain contribu- 
tions from the tail electron population thay may be difficult to evaluate. 
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For example, the tall population contribution to the heat flux may not have a 
simple representation as the gradient of a temperature, a local representa- 
tion, but might be more accurately represented globally (nonlocally ) . The 
reader is referred to the papers of Scudder (1979) and Scudder and Olbert 
(1979) for a discussion of local versus global processes for the non- 
Maxwelliar. solar wind electrons. 

In order to derive a thermal electron heating rate, it is necessary to 
split up the electron population into a thermal (Maxwellian) population with 
energies E < E^ and an energetic (tail) population with energies E > Erp, where 
E T is the transition energy between the two populations. One could derive 
spearate fluid equations for the thermal and energetic electron populations 
with integration limits 0 to Erp for the ambient electron quantities, and Erp 
to 00 for the energetic electron quantities. We adopt this scheme only for the 
thermal electrons treating them as a fluid, and solve kinetic equation 
numerically for the energetic (tall) electron disti, > lution function which is 
to be added to the Maxwellian distribution function. 

The thermal electron energy equation is given by Eq. 2 plus additional 
transport terms on the left side, and with < > modified to be the integral 

over all velocities such that E < E^,. The additional terms are approximately 
given by: 


2E n 




m 


(E - 3/2 kT ) / d ft u f(E T , u) [jf - • u 

Wa A A + jl. A 

+ - u (/ u + ) - ? • u ] , 

3x g m a a m J * 
a e 


(2A) 


where ft is the solid angle of unit vector u, f ( E^, , u) is the electron 

distribution function evaluated at E = E™, and ? is the force per unit 

i m 

mass. The approximation which gives the simple fcrm for the additional terras 
(Eq. 2A) is to equate the density, temperature, and average velocity which 
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characterize the Maxwell distribution function f^, with the same parameters 
evaluated from the full distribution function (f M + f^) for energy E < E^. 

This is a good approximation for most of the ionosphere where kT e « E T , and 
f T « f M for E < E r 

We are interested in this paper primarily in the terms on the right side 
of Eq. 2, the collisional terms which lead to the thermal heating rate. We 
will not discuss further the time and space derivative (transport) terms which 
appear on the left side (Eqs. 2 and 2A). 

With the partitioning of electrons into two populations, we rewrite the 
right side of the energy equation as the sum of two terms, Q e the thermal 

heating rate, and L e the thermal cooling rate. The thermal electron heating 

rate can be written as: 

Q e - / dv 0 (E T -E)(E - | kT e )p + / d v 0 (E t -E)(E - | kT e ) ||* eC (3) 

where <Sf/<Stll is the electron-electron collision operator and 0 is the unit 

ee 

step function, 0(x) = 1 for x > 0, 0(x) » 0 for x < 0. The energy moments of 
the electron-ion and electron-neutral collision operators yield the thermal 
electron cooling rates which have been discussed in the literature (see Schunk 
and Nagy, 1978). The electron heating rate is derived by evaluating the 
integrals in Eq. 3, In order to proceed with this evaluation, we need the 
production rate of electrons, p, the transition energy E^, an explicit expres- 
sion for the electron-electron collision operator, and the total distribution 
function f. The production rate p is obtained during the numerical calcula- 
tion of the energetic electron distribution function. The value chosen for Erj, 
will affect the magnitude of the individual terras in the heat balance 
equation, Q e , L g , and the transport terms. The nominal value for E T is the 
energy where f^ crosses f^ 1 . This value usually has the advantage that: 1) it 
satisfies E T » kT g so the moments of f M over the limited energy range E < E T 


6 



# 


o 


are nearly equal to the moments over the entire energy range; and 2) for E < 
E t , f T « f^ 1 so that the cooling rates are nearly Independent of the tall 
electrons. 

An earlier derivation of the heating rate (Hoegy, 1977) used an electron- 
electron collision operator with simple Debye screening (Shkarofsky, 1961). 

The consequence of that derivation is that it led to an energetic-electron to 
ambient-electron loss rate, dE/dt, having the form of the Butler and 
Buckingham (1962) rate. The more exact energy loss rate of Itikawa and Aono 
(1966) which includes collective wave-particle and wave-wave effects, was 
substituted for the approximate loss rate without being included in an ab 
initio manner. The derivation presented here includes collective plasma 
Interactions ab initio by employing the electron-electron collision operator 
of Kihara and Aono (1963), using their unified theory. Schunk and Hays (1971) 
demonstrated that use of the Butler and Buckingham (1962) energy loss rate 
would lead to as much as a 70 % error in the loss rate calculation when 


compared with the exact loss rate formula containing collective 

interactions. Later, Swartz and Nisbet (1972) asserted that the heating rate 

is not very sensitive to the specific electron-electron loss rate used. 

However, the formula for the heating rate and the calculated distribution 

function both depend on the form chosen for the collision operator. 

6f 

The electron-electron collision operator, is given in Fokker-Planck 
form (Kihara and Aono, 1963) as the first two terms in a velocity-change 


expansion of the Boltzmann collision operator 

5f 3 /oAv v ^ , 1 3 3 . t, Av Av v ^ ,, s 

St ■ '7? • l<4t > f J + 2TJ"; ! f J> <4) 

where the explicit expressions for < —■ > and < — > , are given in Kihara 

and Aono (1963). Their unified theory adds together the impact (binary) and 


collective interaction terms to obtain a non<-divergent electron-electron 
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collision operator. As a result of adding together the binary and collective 
interactions, an exact, convergent term logarithmic In the density is 
obtained. The unified theory has been successfully applied in the calculation 
of transport coefficients (Kihara et al., £963; Kihara, 1964; Itikawa and 
Aono, 1966; and Daybelge, 1969), however, an explicit formula for the 
electron-electron collision operator has not appeared in the literature. 
Appendix A sketches a derivation of the formula for the collision operator, 
using the unified theory. In the next section, this formula is used to derive 
the heating rate expression. 

3. THERMAL ELECTRON HEATING RATE 

The thermal electron heating rate was given in the previous section as 
the sum of two terms. The first term is a direct contribution from electrons 
produced by ionization of neutrals by solar EUV and particle precipitation, 

* J dv e (E t ~ E) (E -|kT e ) p. (5) 

This term needs no further mathematical transformations, and is evaluated by 
numerical Integration of the primary production rate, p. 

The second term, an integral containing the electron-electron collision 
operator, represents heating of thermal electrons by energetic tall electrons, 

q( 2 > = / dv 0 <E t - E) (E - | kT e ) § « ee (6) 

We simplify the form of the second term using the (unified theory) 

collision operator, as given by Eq. A26. Since electron-electron collisions 
conserve energy and mass (these moments of the collision operator over all 
velocity space are zero), we rewrite the collision integral, Eq. 6, as the 
negative of the integral over energies E > E^. and use the collision operator 
evaluated for velocities large compared to the electron thermal velocity: 

Q* 2> - - dir / dv (E - | kT,) { }, (7) 
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where the curly bracket represents the expression in curly brackets in Eq, 


A26. Further integration by parts yields the formt 




S + r m vdv (- f) i&L, (8) 

« V m e (it in 

l 6 

where dE/dt is the test electron energy loss rate given by Eq. A2A and where S 


summarises all the terms evaluated on the energy surface E « E T and is given 


by, 

S«(E T - | kT e )4*Y [l° o f + J - | E T 4irY(I® + J^)f <9A) 

where Y and the integrals 1° and are defined by Eqs. A16-A19. Note that 
in deriving Eq. 8 we have retained only the dominant logarithm terms and have 
neglected derivatives of the velocity-dependent logarithm since they are an 
order of magnitude smaller. The surface term Eq. (9A) is valid for scattering 
particles having an arbitrary distribution function (the argument of the X and 
J integrals). 

Next we evaluate the heating rate, using the collision operator 

given by Eq. A29, valid for Maxwellian scatterers. The scattering particle 

distribution contributes to dE/dt, I, J, and these formulas; the 

scattered particle distribution appears explicitly in the formulas. The 

result has the form of Eq. 8 with the surface term given by: 

E df 

S MT - 4, F 2 (v) [m e - f> £ + < E T - I k V 4 dV !> (n) 

dE ® 

where ^(v) is an integral of given by Eq. A28. The superscript MT 

signifies Maxwellian scatterers and tail population scattered particles. Use 

dE 

of the exact energy loss rate as given by Swartz et al. 1971 in Eqs. 8 and 
9B provides an accurate evaluation of the heating rate for the MT contribu~ 
tion. For the TM and TT contributions (having tail electrons as the 
scatterers), the surface term of Eq. 9 and the ~ expression of Eq. A24 must 


be used 


The thermal electron heating rate is the sum of the three terms given in 
Eqs. 5 and 8-9. The thr^s terms were shown in the review paper by Schunk and 
Nagy (1978) and were first given by Hoegy (1977), however the early form of 
the surface terra was an approximation of the complete form given by Eq. 9. 
Previous expressions for the heating rate have used only the Integral term of 
Eq. 8, and have uBed a variety of lower bounds on the integration. 

It was suggested by Hoegy (1977) that E-p should be the energy at which 
the calculated energetic (tail) distribution function, f T , crosses the 
Maxwellian distribution function, f M . It is not necessary that E T have this 
crossing value, however it is necessary that all cooling terms in L e be 
evaluated over the energies 0 to E.p, i.e. that the .same energy range be used 
for both heating and cooling rates. As E T is varied, both the heating rate Q e 
and cooling rate L fi will change as will the transport terms. If the majority 
of the electron population is the Maxwellian population and if E^ is nearly 
equal to the crossing energy or greater, then the left hand side of the heat 
balance equation, Eq. 2, will be nearly constant and thus the sum Q e + L e will 
be nearly independent of Em, even though the individual terms are dependent on 
E-p. The dependence of Q e on E,p will be demonstrated in the next section. 

4. NUMERICAL EXAMPLE OF HEATING RATE 

We give some numerical examples of the ambient electron heating rate, 

Q e . It is first necessary to solve the kinetic equation for the energetic 
electron distribution function. Using the primary electron production rate, 
shown in Figure 1, calculated for the Venus ionosphere at 60° SZA and 180 km, 
we solve for the photoelectron distribution function fP* 1 such that the total 
distribution is the sum 

f = f M + f ph , (10) 
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where is the Maxwellinn distribution function. Krinberg and Akatova (1978) 

showed that Eq. (10) is a good representation of the total distribution 

function in most of the ionosphere when f ph is calculated according to their 

Eq. (3) and the Maxwell distribution is at temperature T e determined by the 

heat balance equation. TCheir method of computing f ph appears to be identical 

with the method used here when a linear collision operator is employed, A 

thorough discussion of the processes important in the calculation of the 

distribution function is given by Ashihara and Takayanagi (1974). In Figure 2 

we plot the calculated photoelectron flux which is related to f p ^ 1 by, 

Ph 4 2 ph 

«<K) - m f( v)» (U) 

e 

and has the units, if electrons/(cnr sec eV), lhe solid line in Figure 2 
is calculated with the quadratic (in <& p ^) collision operator (Eq. A26); 
the plus signs represent calculated with a linear (in <t p ^) collision 
operator jjivsn by the first term on the right hand side of equation A22, with 
dS/dfc p?fen by the Swartz et al. (1971) formula; and the M's represent <t M for 
Maxwellian flux. The quadratic collision operator terms are non-negliglble in 
the vicinity of the crossing of 4> M and 4» ph and have the effect of enhancing 
the calculated electron flux 4> pfl as was first shown by Krinberg (1973). The 
total electron flux Is the sura of the solid curve and the curve given by the 
M's, and is the flux to be compared with measured electron fluxes. The 
enhancement in 4> p due to the quadratic collision operator may be important in 
such comparisons. Figure 3 is an expansion of the plot of Figure 2 for 
energies 0 to 10 eV, showing clearly the enhancement in the calculated flux 
at E < 1.1 eV. The higher solid line in Figure 3 is the sum of the Maxwelliam 
and the calculated fluxes. In Figure 4 we show the fluxes at an altitude of 
200 km. 

In Figure 3, for 180 km altitude, the Venus ionosphere temperature is 
1280°K and the crossing energy is E T - 1.45 eV; while in Figure 4, for 200 km 
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altitude, the temperature is 2500°K and E T “ 2.63 eV. The Venus parameters 
are from the latest model ionosphere of Brace and Theis, to be published. 

For contrast with the Venus results we show the primary electron produc- 
tion rate calculated for the Earth's ionosphere at an altitude of 200 km and 
solar zenith angle of 42° in Figure 5, and using the reference spectrum of 
Hinteregger (private communication). The calculated flux for this primary 
production rate is shown in Figures 6 and 7. The enhancement in 4> due to 
quadratic terms appears at energies of 2 eV and less. For T fi « 2000°K and N e 
* 1. x 10^, the crossing energy is E T - 1.9 eV. The steepness of the 
Maxwellian flux for the sample cases shown in the figures renders the enhance- 
ment in $> due to the quadratic vs. the linear collision operator insignificant 
in these examples. However, in regions of the ionosphere having high tempera- 
ture T e > 5000° and a higher proportion of tail to Maxwellian electrons, the 
quadratic terms will be important. 

The heating rate terms: Q p , the primary production rate term; Q int , the 
integral term; and S, the surface term are evaluated by substituting the 
calculated and the assumed into Eqs. 5, 8, and 9. These terms are 
shown in Tables 1 through IV along with the total heating rate Q e , 

% “ «p + 1 ta t + S > <12) 

and the nominal heating rate Q nom computed from Q^ nt using from the linear 
collision operator (the plus signs in Figure 2) and using the linear form of 
dE/dt. It is clear from the examples that the individual heating rate terms 
vary strongly with energy E T , and that there is a significant difference 
between the nominal heating rate Q nom computed only from the integral of 
dE/dx over energies greater than Erp, and the net heating rate Q e , which 
includes the primary production term and the surface term. Also, the net 
heating rate is relatively insensitive to the value of E T . 
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5. SUMMARY 


A formula for the ambient electron heating rate, Q e , has been derived 
from the electron kinetic equation. The heating rate depends on the transi- 
tion energy, E^, which defines the upper energy of the "ambient” Maxwellian 
electron population. The cooling rate, L e , also depends on E t j however, the 
sum Q e + L e should be nearly independent of E T . The heating rate is the sum 
of three terms, 


Q + S + Q, 


(13) 


'e ^p '‘int* 

where Q p is the heating from primary photoelectrons, S is a term evaluated on 

the energy boundary E T , and Q in( . is the familiar integral term 

ph 


■ si " dE <- §> F < E >- 

T - dE, ph 


(14) 


- E ; « <-g> *<«. 

ph ph T 

where (E) “ vF (E) is the photoelectron flux, // electrons/(cm^ eVs), and F* 3 * 1 
is the photoelectron energy distribution function, // electrons/ (cm^ eV). 

The value chosen for the transition energy, E^, need not be the energy at 
which the Maxwellian and energetic distributions cross, however it must be 
chosen large enough so that the bulk of the electron population has energies E 
< Erp so that the hydrodynamical parameters.* temperature, density, bulk 
velocity, pressure, heat flux are well represented by integrals over the 
limits 0 to E-p. 

The nominal heating rate calculated in many ionospheric temperature 
studies uses only Eq. 14 and does not generally include terms quadratic 
in either in the calculation of <C> or in the evaluation of the integral. 

In many instances the quadratic terms may not be important, and under some 
conditions, Q p and S may not be signficiant compared to Q£ n (-* It Is important 
though to test the complete formula (Eq, 13) in case these terms are 
significiant. Some studies have used the heating rate of Eq. (14) with a 


lower limit of E T « kT e . Such a heating rate may better approximate Q e than 
Qint with E t * transition energy, however there is no assurance the answer 
will be numerically accurate. 

There is a way of avoiding the problem of the sensitivity of Q e to the 
value chosen for E T . This is to let E T + » and treat the electrons as one 
population. In this case, S + 0 and + 0 so that Q e +■ Q p , which becomes a 

large heating rate; then L e also becomes large. One gains by not having to 
evaluate the surface term and not having to perform the integration in Q£ nt ; 
however the cooling rate calculation requires the addition of integrals using 
the calculated with the Inelastic electron neutral excitation and 
ionization energy loss rates. Hence there is no way of avoiding, first the 
calculation of the energetic tail distribution, and second of performing 
integrals with the distribution, in a study of the electron heat balance 
equation, and in the determination of electron temperatures. 
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Figure Captions 

Figure 1. Primary photoeiectron volume production rate p calculated for the 
Venus ionosphere at 180 km and 60° solar zenith angle. Production 
rate units are // electrons/ (cm 3 eVs). 

Figure 2. Photoeiectron flux 4> calculated for the conditions of Figure 1. 

Plus signs represent the flux calculated with the linear collision 
operator; the M's represent the assumed Maxwell flux; the lower 
solid line is the flux calculated with the quadratic collision 
operator; and the upper solid line is the total electron flux. Flux 
units are // electrons/cm 3 eVs. 

Figure 3. Expansion of the fluxes of Figure 2, for energies 0 to 10 eV. 

Figure 4. Electron fluxes calculated for the Venus ionosphere at an altitude 

, of 200 km and solar zenith angle of 60°. The symbols are the same 

as Figures 2 and 3. 

Figure 5. Primary photoeiectron volume production rate calculated for the 
Earth's ionosphere at 200 km and 42° solar zenith angle. 

Figure 6. Photoeiectron flux calculated with production rate of Fig. 5 and T e 
= 2000°K, N e = 1. x 10 5 e/cm 3 . 

Figure 7. Photoeiectron flux calculated with conditions of Fig. 6, but with N e 
= 1. x 10 6 . 
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TABLE I 


e t 

% 

S 

Qint 

Qnom 

Qe 

1.0 

24 

-454 

3572 

3549 

3142 

1.5 

55 

551 

3264 

3246 

3853 

2.0 

79 

744 

2984 

2966 

3805 

2.5 

109 

908 

2725 

2707 

3741 

3.0 

139 

1038 

2485 

2467 

3662 


Venus calculation at 180 km, 60° solar zenith angle. T^ * 128Q°K, N e = 1.3 x 
10^ e/cm?* Photoe^.ctron and Maxwellian distributions cross at E T » 1.45 
eV. Qp « primary production rate term, S * surface term, Q^ nt » integral term 
using full quadratic collision operator, Q nom " nominal integral term using 
linear collision operator, « net electron heating rate * Q p + S + Q int . 

The units of the heating rates are eV/cm^/sec. 
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TABLE II 


f 


e t 

Q P 

S 

Qint 

dP 

o 

3 

Qe 

2.0 

11 

0 

702 

691 

714 

2.5 

17 

150 

646 

637 

813 

3.0 

23 

197 

595 

587 

815 

3.5 

29 

230 

548 

540 

806 

4.0 

34 

253 

504 

496 

791 


Venus calculation at 200 km, 60° solar zenith angle. T e “ 2500°K, N e = 1. x 
10^ e/cm'*. Photoelectron and Maxwellian distributions cross at E T ■ 2.63 eV. 
Terms defined as in Table I. 
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TABLE III 


e t 

% 

S 

%nt 

^Vtom 

Qe 

1.5 

189 

-297 

1313 

1235 

1204 

2.0 

218 

-134 

907 

878 

991 

2.5 

289 

436 

676 

664 

1400 

3.0 

413 

372 

519 

509 

1304 

3.5 

603 

290 

393 

386 

1286 


Earth calculation at 200 km, 42° solar zenith angle. T e ■ 2000°K, N e » 1.x 
10^ e/crn^. Photoelectron and Maxwellian distributions cross at Erj. ■ 1.9 eV. 
Terras defined as in Table I. 
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TABLE XV 


E>p 

% 

S 

Qint 

^nom 

Qe 

2.0 

218 

791 

2765 

2696 

3774 

2.5 

289 

1226 

2336 

2277 

3852 

3.0 

413 

1166 

1958 

1907 

3537 

3.5 

603 

1123 

1629 

1589 

3354 

4.0 

777 

1184 

1358 

1324 

3318 


Earth calculation at 200 km, 42° solar zenith angle. T e » 2000°K, N e * 1. x 
10^ e/cm'*. Photoelectron and Maxwellian distributions cross at Ey ■ 2.61 eV. 
Terms defined as in Table X. 
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Appendix A 


The explicit form of the electron collision operator is derived using the 
unified theory of Aouo (1962), which unifies the impact and wave collision 
opertors to produce an exact Coulomb logarithm term. In Kihara and Aono 
(1963) Kihara et al. (1963), Itakawa (1963), and Kihara (1964) the exact 
logarithm terra was derived for the energy loss rate of a test electron, the 
diffusion coefficient, and other transport coefficients using the unified 
theory, but the electron collision operator was not evaluated. Shkarofsky et 
al, (1966) developed a general form for the electron collision operator in a 
vector spherical harmonic expansion, but employed approximate logarithm 
terms. The derivation sketched here uses the unified theory and some of the 
methods of Shkarofsky et al. (1966), presents a straight forward method of 
evaluating the impact and wave collision integrals, and obtains formulas used 
in the derivation of the ambient electron heating rate. 

The electron collision operator is taken in Fokker-Planck form as the 
first two terms in a Taylor series expansion of the Boltzmann collision 
operator. The expansion parameter, Av, is the velocity change due to 
collisions. Only terms in the first and second powers of Av contribute to the 
Coulomb logarithm; terms with higher powers of Av are nondivergent and also do 
not contribute to the transport coefficients. The collision operator for 
electron collisions with species s is; 


(*L) 

-'es 


9 v 


(< Av> ^ 
6 £ ) 


3v9v 


(Al) 


where the velocity change — s is given in Kihara and Aono (1963) by their 


wlffi 


Eq. 2.3 for Impact theory, and Eq. 3.25 for wave theory. The tensor term Is 
given by their Eq. 2.2 for the impact theory and Eq. 3.24 for the wave 
theory. Evaluation of the vector and tensor terms is outlined below. 

In the unified theory, the integral over impact parameter, b, has the 
convergence factor exp (-b 2 /2b 0 2 ), where b Q is a parameter obeying the 
inequality, 


close impact radius « b Q « Debye length. (A2) 


These inequalities of b Q with the close impact radius and the Debye length 
allow us to make the mathematical simplifications that lead to the logarithm 
terms. This is illustrated by the Impact parameter integral for the vector 
term: 


2 « f 


bdb e 


j / 2 b 2 /b o 2 


b 2 + , ee s ,2 

' r 


- - it e ? E ± (-5), 


(A3) 


M s 8 


where £ is proportional to the square of the ratio of the close impact radius 


to b 0 : 


o 9 

£ « <_s_) 2 /2b o 2 , (A4) 

V 

Since £ « 1, the exponential integral is well represented by the first term, 
(-£) « In y£> In an expansion in l/£. Thus the Impact part of the vector 
collision term is: 


At impact 


g 


Y5* 


(A5) 
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where K ■ 4n e 2 e 2 /m 2 , g * v - v is the relative velocity, and y is Euler's 
s e s 

constant, y * 1.781072. 

Integration of the tensor term also Involves th? exponential Integral and 
yields: 


<AvAv> g 
At impact 


K / d \ {(lt " 1 ln 7? 

+ (3 gg - It) I), 


(A6) 


where tt is the unit tensor, and g » g/g. 

The wave (collective interaction) integrals are more complicated than the 
impact integrals. We consider first the tensor term (Eq. 3.24 of Kihara and 
Aono, 1963) and use the order of integration of Shkarofsky et al. (1966). The 
Dirac 6 function in the to integral sets co » it • v and simplifies the 
dependence of the integrand on the magnitude of it; with the unified theory 
convergence factor, exp (-k^b^/2), the k integral gives: 


<AvAv> 


At wave 2 tt 


‘hi! 


dv f (v ) 
s s s 

8 


A A AAA 

dk 6 (k • g) kk 


[ (cos y - ~ sin y ) ln ( tt — z ) 

y y/x 2 +y 

+ (sin y - ~ cos y) arc tan (■£)], 


where, 


b 2 k 2 
x = J - -- -- - R (x g ) 
s 

b V 

y - l — t - 1 < x *> 


(A7) 


f 2 t- 2 

R(t) = 1 - 2t e" r ' f dz e Z 
_ t 2 

I(t) = /iTt e __ 

x s - * • y 'fzk 
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2 


2 4 * "s e B 


kT 


9 9 

and b Q k g << 1 from the inequality, Eq. A2, where k g is an inverse Debye 
length. From this inequality, x « 1 and y << 1, thus the square bracket 
expression in Eq. A7 reduces to 


In \ r - ~ arc tan (•£) , 

Y /x 2 + y 2 y x 


the same functional form obtained by Daybelge (1969). The remaining integral 

A A A 

over all angles of k such that k • g = 0 is similar to the integral over 

A A A 

impact angles b such that b • g = 0, however, the wave case is more 

A A 

complicated because of the dependence of x and y on k. The k integral has the 
form: 


A A A 


/ dk <S (k • g) k k [ ] 


= ^ 2lTd 4> ((?t - gg) cos 2 (j> + yy (1 - 2 cos 2 .)))) [cos 2 <|)J, 


(A8) 


where x, y, and g form a right handed orthogonal coordinate system with y in 

the direction of v x v , cos^ = k • x, and [cos .)>] indicates that the square 

s 

bracket expression is an even function of cos^. 

We are interested in evaluating the collision operator for energetic 



electrons, when v \J — » 1; in this case [ ] has a logarithmic dependence 

2 s 

2 S 
on cos i)>, and the k integral reduces to, 


+ + A A 1 a * A 

2v (11 - gg) ~ In j - 2 7T yy, 
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2 

where A represents the argument of the logarithm multiplying cos <{>. Therefore 
in the argument of the logarithm, cosif) is replaced by — using the integrals: 

0 /^ n d<}> In (cos%) ■ 2tt In (— ) 

o / 2tt d<J> In (cos 2 ^) cos 2 (J> - 2ir ^ (In (-|) + 1). 


4- + AA 

Eventually, in the integral over v with fixed v, the yy term becomes 

s 

1 AA 

Y (11 - vv). Therefore the tensor wave integral is: 


<AvAv> 


s 


At wave 


» K / dv 


? fsllal ( tr^l 

s g 2 yn 2 


], (A9) 


where, 


.2 2 

b w 

n * _2_ _E 

n 2 2 1* 

V 1 4 


o n e 

2 r / s s 
w = ) 4u 

P s 


m 


vv 


5- /1-z 2 , 


Z * V • V 

s 

The evaluation of the wave vector term follows the above steps, except 

that the derivatives — and (from Eq. 3.25 of Kihara and Aono (1963)) must 

3v 3v 

be transformed using the relation 


lc — — 6 (u>-* & • v) *k 2 —6 (& » v — w) 

3v 3v 
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and integration by parts with v . The result is; 
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<Av> ra +m . a , . . , 

-TFwave ' - K "V / *.*.<*.> 7 I IT Tn - ■ (A10) 

8 3g 


The derivative, — , acting on q produces nonlogarithm terras which we neglect, 

Sg j 

therefore the derivative acts only on — . 

g 

Summing the impact and wave contributions gives a net logarithm argument 


which is independent of b Q : 


1 . 1 , 1 . 1 , . 

■x In — + t In — = In A, 

2 Y5 2 yn * 


where 


Y ee ii) 

' s p 


(All) 


Thus the unified theory vector and tensor collislonal velocity changes are; 


<Av> m -fm + 

-ST 1 = K -V s ! 4 Vs ( *s ) % [In A - VjJ , 

S g 

<AvAv> f s (v ) a. .. ^ , 

— ££ = K / dv g — [(ft-gg) In A - (tt-vv) 1 ^ 


( A12) 


+ (3 gg - tty/rf. 


(A13) 


This result, Eqs. A12 and A13, is valid when v is much greater than the 

thermal velocity. This restriction can be relaxed, however, the more general 

form is not needed here. In the following, we further develop the vector and 

tensor terms using the techniques of Shkarofsky et al. (1966). 

The integration over v is carried out by expanding f (v ) and the 

s s s 

integrand in vector spherical harmonics. A convenient means of evaluating the 
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expansion Is to employ the Rosenbluth potentials (Shkarofsky et al. , 1966). 

f* * 

To simplify the derivation we neglect the dependence of A on v • v and assume 

s 

f (v ), the distribution of scatterers is isotropic. Henceforth we retain 

8 S 

only the dominant logarithm terras. 

The results are: 


Y At 


. <AvAv> 
_l_ s_ 

Y At 


+• m +m 
v e s o 

3 m o 9 
v s 


It- 


vv ,^.o 


(r + 


( Al 4) 


+ (- 3 ^ v - t — ) <I° 2 + J_°>. 


(A15) 


(4Tr) 2 e 2 e 2 

Y - j—5* ln V 

m 


( A16) 


2„ s v- 


o yee w * 
s p 


. A Vs. 
ln <T ~Z) 


tO 1 fv2j c o, N n , 2 v \ 

1= — f v dv f (v ) v (1+ — z 7 J 

n n A s s s s s v In A ' 


(A17) 


( A18) 


T° - JL 
n n 
v 


i A 2 / 2 \ 

9 „ InC-r v h ) 

v dv f (v ) v (1 + — • » — t ). 

s s s s s InA 


(A19) 


The collision operator, Eq. Al, can now be evaluated using Eqs. A14 and 
A15, and the following transformations: 


<Av> 


7* ( " 

8v 


s . v 1 d , + 

— f e> “ — dv' (v v 
V 


<Av> 

At e'* 


(A20) 


a 2 <AvAv> . . 

; r l f ) * - 1_ ±_ 

^ At e ' 2 dv 

3v3v v 


(v 


< Av • A v> 


At 


f J 


, A , ^ . . <AvAv> 

+ ( — — [v vv : — — f ] 1 . 

2 dv W dv At e ' 

v 


(A21) 


Combining the terms, we find: 
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6f . ) 

6t es 2 dv '■m dt e' 


v 2 dv W dv 


<A?AV > 8 
(v vv: — 



(A22) 


where — ■ is the energy loss rate of a test electron due to collisions with 
particles of specie s, 


dE r + 

dt = m e tv 


<Av> <Av*Av> 

5. + 1 5.1 

At + 2 At J> 


(A23) 


Evaluation of the energy loss rate using Eqs. A1 4 and A15 gives: 


ra dt 
e 


- V (- 


m 

-S. 1° + J 
m o -1 1 
s 


(A24) 


This expression is a generalization of the fast test particle energy loss rate 

of Kihara and Aono (1963) since the full dependence of the integral on v, v g , 

and v • v has been considered in the derivation. The cited work used the 
s 

approximation of replacing g by v. The second term in the collision operator, 
Eq. A22, is: 


1 

2 


v 


vv: 


<AvAv> 
s 

At 


Y 2 
3 V 




). 


(A25) 


thus the electron - s particle collision operator is given by: 


K) 

St'es 


1 


m 


1° Y f 
2 . 1 in 1 o 

v dv s 


+ I Y < l0 2 + 


df 



dv N 


(A26) 


31 


This result is similar in form to Eq. 7—7 lb of Shkarofsky et al. (1966), 
differing only in that the present form contains the exact Coulomb logarithm. 

In the course of the derivation we have made simplifications which amount 
to neglecting the non-logarithmic terms and we have assumed v » the electron 
thermal velocity. These simplifications can be relaxed at any stage in the 
above derivation; also the case of anisotropic distribution functions can be 
included. 

The argument of the logarithm in the final results Eqs. A16 and A17, is 
identical with the logarithm of Kihara and Aono (1963), but differs from the 
results of the later paper of Itikawa and Aono (1966) where arbitrary velocity 
v was allowed. In the latter, the logarithm is split into a term independent 
of v and a complicated function, G, containing the v dependence. At 
sufficiently high energies, a quantum mechanical evaluation of the collison 
operator is necesary. The quantum mechanical calculation of Kihara (1964) 
gives the logarithm argument: 


mv^/hw . 

P 

The tedious evaluation required of the G function, prompted Swartz et al. 
(1971) to make an analytic fit to the energy lorfs rate ^ for both the 
classical and quantum mechanical formulas, however, their analytic function 
best represents the classical formula with the logarithm argument of Eq. A17 
at all energies and does not seem to have been fit tq the quantum mechanical 
result. 

It is useful to derive an alternate formula for the part of the electron- 
electron collision operator having Maxwellian scatterers. We begin with Eq. 
A22, with s = e, f g = Maxwellian distribution, and dE/dt given by the general 
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formula of Itikawa and Aono (1966) or the analytic formula of Swartz et al. 
(1971). Let 


■V' ifr-rwr" 1 


F 2 (v) 


V vv 


+• . 


l_ <AvAv> 
2 At 


(A27) 


A simple equation is obtained for F£(v) from A22 when f e + f fi M (the Maxwell 
distribution function) * 


2 dv 
v 


v dE 

m dt 
e 


£ M + if 

e v dv 


< F 2 f e M)1 


0 , 


since F£(v) * s independent of the form of f e< 
The solution is given by: 


P 2 (v) = 


- / 

v 


'2 ' 
v dv 


m 


dE 

dt 


2kT 


i ' 2 2 \ 

(v - v ) 


(A28) 


Thus the electron-electron collision operator for Maxwellian scatterers is 
given by: 


*L) 

6t 'ee 


1 _ 

2 dv 
v 



m , df 

<v) + 7d^r 1} , 

e 


(A29) 


and this formula is valid for arbitrary velocity v, since it contains the 
exact energy loss rate dE/dt. For tail distribution scatterers, Eq. A26 
should be used. 
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